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Abstract 

We compute Greeks for stochastic volatility models driven by Brownian informations. We use the 
Malliavin method introduced by !1 for deterministic volatility models. 
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1 INTRODUCTION 

The application of the Malliavin calculus to the computations of price sensitivities were introduced by 
PQ for models with deterministic volatility. In this work, we compute the Greeks, for stochastic volatility 
models where the underlying asset price is driven by Brownian information. We consider stochastic 
volatility models, since these models, unlike those with deterministic volatility take into account the 
smile effect. 

Let (B t )t£\o t T} an( A (B t )t£[o.T] De two independent Brownian motions. We work in a filtered probability 
space (fl,^F, (^Ft)telo,T]jP)> where (Ft)te[o,T] is the naturel filtration generated by B and B . We consider 
a market with two assets, a riskless one with price (e-k rsds )tG[o,T]i where r is deterministic and denotes 
the interest rate. And a risky asset (<S't)tG[o,T] to which is related an option with payoff /(St)- {St)te[o,T] 
has a stochastic volatility and is given by 

dS 

— - = ntdt + a(t, Y t )dB t , where alY t = pj alt + a\ ' [pdB t + dB t ], t E [0, T], 
St 

with So — x > and Yq = y G M. \i, fi Y , a Y are deterministic functions, p € R and a £ C 2 ([0,T] x M) 
such that for any t G [0, T], o~(t, .) ^ 0. The market considered here is incomplete. There is an infinity of 
E.M.M -Equivalent Martingale Measure- (i.e a probability equivalent to P under which the actualized price 
(Ste~fo r< * <is )t e [o J T] is a martingale). Let Q be a fixed P-E.M.M. Q is identified by its Radon-Nikodym 
density w.r.t P, denoted Pt and given by 

PT^expN a s dB s + (3 s dB' s -^ (a 2 s + ff)dsj , 
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where (at)te[Q,T] and (A)ie[o,T] are two predictable processes s.t. a t — — ^ryn and P is arbitrary. Now 
let, for any t S [0, T], W t — B t — J Q a s ds and — B t — J* f3 s ds then by the Girsanov theorem VF and 
W are two Q-Brownian motions. In the following we work with a fixed P-E.M.M. Q and we will use E[.j 
(instead of -Eq[-]) as the expectation under the probability Q. We have, under Q, for any t £ [0, T] 

^1 = rt dt + a(t, Y t )dW t , dY t = (tf + ^^f^T + dt + P°Y dW t + ajdw' t . (1.1) 

St V °\P,Yt) J 

2 MALLIAVIN DERIVATIVE ON WIENER SPACE 

In this section we give an introduction to the malliavin derivative in Wiener space and to its adjoint 
the Skorohod integral. We refer for example to [2]. Let {DY)te[o,T\ be the Malliavin derivative on the 
direction of W . We denote by P the set of random variables F : fi — > M, such that F has the representation 



jf h(t)dW U ...,j o fn(t)dw}j , 



F(uj) = fU h(t)dW u 

where f{x±, . . . , x n ) — ^2 a a a x a is a polynomial in n variables x\,...,x n and deterministic functions 
fi £ L 2 ([0,T}) . Let ||.||i, 2 be the norm 

11*111,2 := \\F\\ L * m + \\D w F\\ L 2 { ^ T]xn) , F € L 2 (n). 

Thus the domaine of the operator D w , Dom (D w ), coincide with P w.r.t the norm ||.||i,2- The next 
proposition will be useful. 

Proposition 1 Given F = f (j* fi(t)dW t , ■ ■ ■ , ^ f n {t)dW^ £ P. We have 



fn(t)dW t f k (t). 



fc=0 

To calculate the Mallaivin derivative for Ito integral, we will use the following Proposition. 

Proposition 2 Let {v,t)t£[o.T] be a T \ — adapted process, such that Ut € Dom (D w ), we have 

Df f u r dW r = [ (Df r u 8 )dW s +u t . 
Jo Jt 

From now on, for any stochastic process u and for F £ Dom (D w ) such that u.D w F £ L 2 ([0,T]) we let 

rT 
/o 

Let 6 W be the Skorohod integral in Wiener space. We have 6 W is the adjoint of D w as showing in the 
next proposition, moreover its an extension of the Ito integral 



D?F := (D w F,u) LH[Q , T]) := f u t D? Fdt. 

Jo 
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Proposition 3 a) Let u e Dom (S w ) and F e Dom {D w ), we have E[D™ ' F] < C(u)\\F\\ 1<2 , and 
E[FS w (u)} =E[D^F]. 

b) Consider a L 2 (Q x [0,T]) -adapted stochastic process u — {ut)te[o,T]- We have S w (u) — u t dW t - 

c) Let F G Dom (D w ) and u G Dom (8 W ) such that uF G Dom (8 W ) thus 5 W (uF) = F6 w (u) - D^F. 

3 COMPUTATIONS OF THE GREEKS 

The computations of Greeks by Malliavin approach rest on a known integration by parts formula -cf. pQ- 
given in the following proposition. 

Proposition 4 Let L be an open interval of R. Let (F^)^^i and (H^)^i, be two families of random 
functionals, continuously differentiable in Dom (D w ) in the parameter £ G /. Let (ut)te[o,T] be a process 
satisfying 

D^F^ ? 0, a.s. on {d c F< ^ 0}, C G /, 
and such that uH 1 ^ d<^F^ / F^ is continuous in C, in Dom (6 ). We have 

±E [H<f (2*)] = E [f (FC) (^^(u) - D? {^0-) + 
for any function f such that f G L 2 (f2), £ G J. 



Our aim is to compute the Greeks for options with payoff /(St), where (St)te[o,T] denotes the underlying 
asset price given by We have 

Sr=xexp(j o a(s,Y s )dW s + (r s -^a 2 (s,Y s )ds)j , (3.1) 

Let £ be a parameter taking the values: the initial asset price x — So, the volatility a, or the interest rate 
r. Let C = E[f(Sf)} be the price of the option. We will compute the following Greeks: 

dC d 2 C dC dC 

Delta = — — , Gamma = -?r~o"> Rh° = tt~ > Vega = — — , 
Ox Ox z or oa 

The next Proposition gives the first, second and third order derivatives of St w.r.t D , needed for the 

computations of the different Greeks. 

Proposition 5 For < t < T, we let G(t, T) := a(t,Y t ) + jf ^(v,Y v )DfY v (dW v - a(v,Y v )dv). We 
have DY St = STG(t,T) and thus 

D^S T = S T / u t G(t,T)dt (3.2) 
Jo 

D^D^St = s t((J u t G(t,T)dt\ + J J u s u t D s G{t,T)dtds\ (3.3) 
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n W n W r,W o 

D u U u D u b T 



u t G{t,T)dt 



[ [ u s u t D s G(t,T)dtds 

JO J s 



u t G(t,T)dt 



+ 



r T r T r T \ 

I]] u r u s utD r D s G(t,T)dtdsdr 



where 



DYG(t,T) = ^(t,Y t )D?Y t + [ T DYY v D?Y v 
dy J t 



d_ 

dy 



da 
dy 



(v,Y v ){dW v ~a( Vl Y v )dv) 



D^DfG{t, T) 



+ f ir( v > Y v )DY Df / Y v (dW v - a(v, Y v )dv) 
Jt oy 

£ ({d^dYy v dYy v + dYy v d^dYy v )^{v,y v ) 
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+D? Y V D™ Y v D™ Y v — (v,Y v ) 

+^(v,Y v )DYDYDfY v + d ^{v,Y v )DfY v DYDfY v ) dW v 



{DYDYY V DYY V + D^Y V D^D^Y V )^ ( a(v,Y v )^(v,Y v ) 



lW V r,W nff 







dy 



da 



dy 



±(a(v,Y v Av,Y v )\ /),"}:;/)!' /)," Y, 



dy 



dy 



da 



dv 



a(v, Y v )^(v, Y V )DYDYDYY V 
dy 

+ d ^Y t )DYY t DfY t + ^.(t,Y t )D?D7Y t . 
Proof. By the chain rule of Dj and thanks to Proposition [2] we obtain 



DrS T = S T Dl 



a 2 (.s,n) 



ds + Df a(s,Y s )dW^J 



£ D w ^ 2 {^Ys) ds + £ D w a{Sj Ys)dWs + ff(tj Yt) ^ = StG{ ^ t)) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



which gives (13. 2|) . (|3.3p and (|3.4p are immediate by the chain rule of D w . Concerning (|3.5p and (|3.7p we 
have for < r < s < t < T 

DfG{t,T) = ^( t ,Y t )DYY t + £dY (^(v,Y v )DYY v y Wv 

- [ Df (a{v,Y v )^{v,Y v )DfY v \ dv 

= ^(t,Y t )DYY t +£ (dYy v dYy v ^{v,y v ) + ^{v,Y v )DYDYY v ^ dW v 



DYY v DfY v ^ ( a (v,Y v )^-(v,Y v ) ) + a(v,Y v )^(v,Y v )D™ D?Y V 



d 



dy 



da 



dy 



da 



dy 



dv. 
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And 



D?DYG{t,T) = j\^{DYY v DYY v d ^Y v ) + ^-{v,Y v )DYDYY v 



I 



D 



w 



d 



dy 



da 



D^Y V D^Y V — a(v,Y v ) — (v,Y v ) 



dy 



da 



+a(v, Y v )^(v, Y v )D^DfY v 



dv 



2 a 



+^(t,Y t )D^Y t DYY t + ^(t,Y l )D) ] Z?» V, 



{D™DYY V DYY V + DfY v D™ DYY v y—{v, Y v ) 



dy 



+D? Y v Dg Y v Dt Y v ^(v, Y v ) 
da, , r w w w , r d 2 - 



{v,Y v )D™D™D?Y v + ^—(v,Y v )D^Y v D^D^Y v ) dW v 



dy 



dy 



d 



{D™DYY v DYY v + DYY v DfDfY v )^- a(v, Y v )^(v, Y v ) 



dy 



da. 



dy 



+ 



dy 2 
d 



a(v,Y v )^-(v,Y v )) D)'Y r Dl ] D," V, 



dy 



da 



dy 



( a{v,Y v y—{v,Y v ) I DfY v D™DfY v 



da 



dv 



+ a(v,Y v )^-(v,Y v )D^D^D^Y v 
+ ^(t,Y t ) D w Yt DYY t + 9 ^{t,Y t )DYDfY t . 



Concerning D^Y V , it can be explicitly computed when (3 V = (3(v, Y v ). We have for < t < v < T 



D W Y - D w 
u t i v — u t 



f 

JO 



Y Y ' ol fa 
Ma + ° a 



a(a,Y a ) 



+ (3(a,Y a )a Y ^j da + pa Y a dW a + a 



a dW a 



= £ D? (/£ + al + fta, Y a )at) da + D? jT P a Y a dW a 

= [ olDf (^=^ + /3(a, Y a )j da + jT D? pa Y dW a + pa Y t 
r a -p a df3{a,Y a y 



□ 



D W Y da 

) t "° \a 2 (^Y a ) dy J Vt y « da 

So for t fixed in [0,T], the Malliavin derivative of Y v for v E [t,T] : (D^yY v ) ve [ t T ^, satisfies a stochastic 
differential equation, its solution is precisely 

n ff v Y ( [ V ( Y r a-V>a d0(a,Y a )\ \ 

D t Y v = pa t exp J- ^ ^ « € [i,T]. 

In this case for < r < s < t < v < T, Df DfY v and Df DfY v are as follow 



D w D w Yy = BfY v D 



w 



r a -p a d(3(a,Y a ) 



a 2 (a,Y a ) 



dy 



da 
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D w D w D w Ys 



D?Y, 



' " < t \r" a a*(a,Y a ) • dy 

r a -(J. a d 2 f3(a,Y a )\ w 



{D™DYY v ) I ( 2*1 



t \ 2aa ^{a,Y a ) ' 8y 2 



D^Y a da 



Y r a Ha d 2 /3( a ,Y a )\ Ty 



oH a ,Y a ) dy 2 



D™Y a dv a 



-D W Y 

u t i v 



+DfY v 



6(7. 



2a 



y r " ~ Ma _ d 3 l3{a,Y a ) \ w w 
a 4 (a,Y a ) dy 3 > r a s 

r a -fi a d 2 f3(a,Y a )\ w w 



a3(a,Y a ) dy 2 



D™D™Y a da. 



3.1 Delta, Rho, Vega 

The Delta, Rho, Vega can be computed using a first order derivative of C = E[f(S T )] w.r.t (. We have 
using Proposition 2] 



d_ 



E 



/(4) 



= E 



-s w {u)-d: 

T 



'u 



Next we compute the Delta, the Rho and Vega can be computed by the same way. The Delta corresponds 
to C = x, so 8qSt — d x Sr — jSt and we have 



Delta 



E 

E 
1 



f(Sr) 



f(Sr) 



(A 



S- 



St 



5 w {u)-D 



w 



8 w {u)-D\ 



St9 x St 

Srp 



E 



f(S T ) 



5 w (u)~2S T + 



Jo u t G(t,T)dt 



S T D™D?S T 



= -E 

x 



f(S T ) 

\£u t G(t,T)dt 



S w (u)-2S T 



S T Jq utGjt, T)dt) 2 + St Jq u s D s (J T u t G(t, T)dt)ds 



J* utG(t,T)dty 



1 



-E 



f(Sr) 



S w (u)-S T 



Jo u t G(t,T)dt 



' x _ J T UsjJ u t D s G{t,T)dt)ds 
J T u t G(t,T)dt) 2 



where DfG{t,T) is given by (JH 
3.2 Gamma 



The Gamma is computed using the second order derivative of C = E[/(St)] w.r.t x given by 



^E If (S T )} = |-Delta = -%-E [f(S T )H] 

OX OX X ox 



1 

—E 

x 



ft a \ ( Hd xS T f. Wf ^ w f Hd xS T 



d x H 
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= l E 

x 



HS T , w , s (H{D^S t ) 2 + StD^StD^H^HS t D^D^St 



+d x H)] , 



where 



/o u t G(t,T)dt \ (/„ u t G{t,T)dt) 2 

w„ £u s {fu t D s G{t,T)dt)ds xW( 

u tl = 7p (Ml 

u t G{t,T)dtf 

fo u s (fi u t DsG{t,T)dt)d. 



St ( I U t G(t, T)dt 

\Jo Jo u t G{t,T)dt 

/ T Ur / r T Us (/J u t D r D s G(t, T)dt)ds (lo u *(£ u t D s G(t, T)dt)ds 



(J T u t G(t,T)dt)2 ' (Jo u t G(t,T)dty 

l e (, u s (£ Ut D s G(t,T)dt)ds' 

T 



d x H = --S T I 1 - — v . 

* V (Jo u t G(t,T)dt) 
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DfG(t,T) and DfDfG{t,T) are given by formulas ([23]) and (1X71) . 
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